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SUMMARY AND INTRODUCTION 



This report presents the results of experimentation with a single 
jet-augmented flapped airfoil for generating gusts of high reduced 
frequency for the range 



in a wind tunnel. The experimental data from a single jet flap in 
conjunction with the theoretical assumption of superposition, are 
used to postulate combinations of trailing vortices from an assumed 
array of multiple jet flaps for the production of an oscillating 
cross -sectionally-uniform gust. The objective of this report is to 
study and understand the effects of different parameters such as flap 
deflection angle, flap frequency, dynamic pressure, and the position 
at which the induced velocity is measured, on the uniformity of the 
generated turbulence across the wind tunnel test section and finally 
to combine these system variables to produce sufficiently uniform gusts. 

Three areas of work were the main points of concentration: 

1. ) Experimental testing in both the static and dynamic cases 

with a single jet flap or vortex generator in the wind 
tunnel. 

2. ) Superimposing the results of the experimental data into an 

array of four generators in an attempt to produce a uniform 
cross-sectional gust. 

3. ) Comparison of the gust model designed using the experimental 

data and superposition with the model using a theoretical 
analysis of a single trailing vortex and superposition. 
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The testing was accomplished in the Princeton University 4 1 X 5 f 
wind tunnel with jet augmentation of a simple plain flapped airfoil. 

The jet augmentation increased greatly the magnitude of induced velocity 
generated, and the use of a flapped airfoil decreased the size and inertia 
of the moving components. 

The experimental results consist of measured values of induced 
velocity due to a single half span vortex generator. Both static and 
dynamic data were taken. In the static case, the magnitude of induced 
velocity was dependent on certain variables such as flap deflection angle 
sensor vane position in the test section and free stream velocity. 
Statically, many combinations of parameters were observed and recorded. 
Thus, the partial derivative of induced velocity with respect to each 
parameter was separately observed. For the dynamic data, induced 
velocity was dependent on- Strouhal Number. Bi-directional as well as uni- 
directional gusts were generated and were compared. An analysis of the 
effects of high frequency cancellation and vortex filament movement 
illustrated a similarity in peak to peak magnitude in both types of 
generated turbulence. 

The second area of study was the use of the experimental data in the 
hypothesis of an array of generators that would produce the desired uniform 
gust model. Superposition was used to pattern the symmetrical array of 
trailing vortex filaments into this desired model. 

Finally, a similar array was developed from a trailing vortex 
filament in an inviscid flow field. This theoretical model was compared 
to the model using the experimental data and superposition previously 
designed. They were in close agreement. 
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DESCRIPTION 



Gust simulation in the 4X5 foot Princeton University wind tunnel 
was conducted using a blown flap gust generator. No attempt was made to 
improve the design of the gust generator based on experimental observations. 
However, changes were made in the sensors during the trial runs to improve 
the data output and to increase experimental accuracy. The individual parts 
will be discussed separately including the changes in the sensor construction 
during this trial run period. 

Wind Tunnel 



The tunnel test section (Figure l), nominally four feet by five feet, 
is forty-two inches high and fifty-nine inches wide. This Princeton 
facility has the least tunnel turbulence of all the wind tunnels at the 
Forrestal Campus. Its dynamic pressure is variable over a range of tunnel 
velocities up to approximately 100 mph. The maximum tunnel free stream 
velocity used in the experiments was 100 ft/sec which is below the upper 
limit where tunnel turbulence tends to be amplified. The test section 
has a streamwise length of 80 inches, thus allowing surveys at three positions, 
two, six and ten chord lengths downstream from the jet flap model. 

Inherent in any wind tunnel is its turbulence caused by wall surface 
irregularities and friction (shear) stresses. This turbulence has its own 
characteristics of frequency and amplitude which must be filtered out of 
the final data in an attempt at some degree of accuracy. The frequency 
of the turbulence was recorded in the initial trial runs as 40-50 Hz, 
which was well above the testing range of the experiment and was easily 
distinguishable and filtered. The magnitude of the turbulence varied 
between one and two degrees of free stream deflection in a steady state 
condition, and wa.s filtered out of the data presented. 
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JET FLAP VORTEX GENERATION 



The jet flap was designed specifically to increase the lift 
capability of a particular airfoil section. Since thesis increased, 
the circulation also increases and hence the strength of the trailing 
vortex far exceeds that of a simply flapped airfoil. 

To create a sinusoidal gust the airfoil section was chosen as a 
NACA 0012 symmetric airfoil. The flap, (Figure 2) was taken to be a 
,45c flap for reasons discussed in APPENDIX I. Blowing over both sides 
of the flap was necessary to insure a symmetrical gust profile. 

This trailing vortex generator was a half span model mounted to 
the tunnel floor plate to facilitate oscillating the flap. Mechanically, 
the flap is oscillated by means of a connecting rod attached to a driven 
fly wheel as in Figure 3- The wheel has different radial attachments 
for the connecting rod. By increasing the radial position the actual 
flap deflection can be increased both positively and negatively. 

The drive wheel is powered by an electric motor through a 10 to 1 
ratio reduction gear box with a right angle power transmission. The 
motor was rated for maximum power output for 27 volts with 69 amperes of 
current. Throughout the entire range of tests the voltage never increased 
over 27 volts with 10 amperes; therefore, the maximum power output was 
never approached. All powered components are controlled from a central 
console as sham in Figure 4. 

Blowing over the flap was- accomplished by the introduction of a' 

.015 in ff width slot over the full twelve inch semi-span of the airfoil. A 
large plenum chamber was designed into the flap construction to feed air 
to this slot. The Schramm Compressor was able to produce a plenum chamber 
pressure of approximately 30 psi, enough to choke the slot and give spanwise 
uniformity to the jet sheet. Line loss was almost 50 psi from the output 
of the compressor to the point of induction into the plenum chamber. 
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FIGURE I. DOWNSTREAM VIEW OF TEST SECTION. 
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FIGURE 2. JET AUGMENTED FLAPj VORTEX GENERATOR. 







FIGURE 3. DRIVE MECHANISM FOR JET FLAP; ELECTRIC D.C. MOTOR, GEAR BOX, 
DRIVING WHEEL, AND CONNECTING ARM. 
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FIGURE 4. CONTROL CONSOLE AND WIND TUNNEL. 
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FIGURE 5. GUST SENSOR 




3 



INSTRUMENTATION 



Air stream following vanes mounted on low inertia potentiometers 
in Figure 5 were used to measure the steady state and time, varying 
response of the free stream air to the sinusoidal oscillation of the jet 
flap airfoil upstream. The dynamics of the vanes and ,their design 
characteristics are discussed in Appendix IV. Adequate speed of response 
and no overshoot in time response are two desirable characteristics of the 
vanes. After many different types of vanes were designed and tested, the 
vanes shown in Figure 5 were used. They have low aspect ratio conical 
bodies with laterally-protruding wings at leading and trailing edges. 

By inducing an impulse perpendicular to the free stream in a steady 
state condition it was observed that with an aft moment arm of 3 inches 
and a forward arm of 1.5 inches, the vanes are critically damped at all 
tunnel velocities. This is advantageous for accuracy in data recording. 

The natural frequency was in the vicinity of 30-35 Hz which is well beyond 
the range of experimentally-induced frequencies. 

The probe shows the configuration of the potentiometers 
and vanes in the free stream. One potentiometer is mounted laterally 
placing its vane in position to record any vertical component of 
deflection in the stream. The vertically mounted pot and vane will record 
the horizontal component of disturbance. By vectorially adding these 
components we arrive at the velocity vector induced at that point. If we 
assume that this velocity vector is perpendicular to a radius to the 
vortex core then this represents a good method for experimentally locating 
the vortex core position laterally in the tunnel. 

Electrically excited by a 15 volts source the potentiometers were 
attached to a Sanborn heated stylus recorder and a Brush pen recorder. 

With maximum frequency capability of 60 cps and 100 cps, the Sanborn and 
Brush recorders, respectively, could easily accommodate even the high 
tunnel turbulence frequencies. Calibration was accomplished by mechanically 
varying the vanes through a known angular distance. Brush and Sanborn 
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recorders had 1 ram displacement on the recorders equal to 1.5 degrees 
of vane deflection for the former and 1.3 degrees for the latter. 

EXPERIMENTAL PROCEDURES 



The experimental testing was accomplished over a period of six weeks. 
A number of assumptions had to be made to simplify reduction of the data. 
Incompressible flow and constant temperature were assumed. With a 
number of variable parameters , the most efficient means of surveying 
cross sectional areas of the 'tunnel was devised. 

The variable parameters consist of the following: 

1. ) Tunnel dynamic pressure which alone varied the value. 

\ C"> ^ 

2. ) Blown flap deflection, 

3. ) Frequency of flap deflection,/,^. 

4 . ) Position at which the induced velocity is measured in 

test section, • 

The technique consisted of the following steps: 

First fix the following: 

r~ 

1. ) The specific -y value. 

2. ) A specific vane position in the tunnel, . 

3. ) The initial tunnel dynamic pressure, q. 

4 . ) The initial flap frequency, 

The range' of these parameters is very important in gust simulation. 
The various combinations can change the magnitude, direction, and span- 
wise distortion of particular gusts. 
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For the present discussion the following magnitudes of the 
parameters were employed. 

The initial condition consisted of 

a) 

(2) VoO = 100 5t/&c 

(3) CC>f O 

First, a static case was run for each new and Sf . The blown flap 
was manually cranked, slowly, through a single cycle. This created 
a good static trace. Then the flap was oscillated through a range of 
1 hz 9 2 hz 9 4 -hz , 8 hz , and 16 hz , with the run time of 10 to 15 
seconds. Next the V£ 0 was lowered to 75 ft/sec and the static and dynamic 
cases were run again. Finally, the ^ was fixed at 50 ft/sec and the process 
was repeated. The V\> and changes were made for all £^- t in a specific cross 
section of tunnel at x = 2.0c downstream from the model. After this cross 
section was surveyed, the vanes were moved to the x = 6.0c downstream 
section and finally, the x = 10.0c downstream section. The jf magnitude 
was then changed to accommodate a variation of ± 45° . After the entire 
tunnel was surveyed again at df = 1 45° , the 5f-was changed to ± o0° 
for the third and final run through the tunnel. During the = ± 30° 
survey, approximately eleven to thirteen evenly spaced points were surveyed 
in each cross section. For the (if = ± 45° and ± 6o° , only four of the 
centrally located points in each section were observed. This reduction in 
survey points does not decrease the accuracy of the experiments. Induced 
velocity magnitude variations are still discernible and comparisons with 
different parameter regimes can still be made. 

> 

A few runs were made without blowing augmentation for comparison with 
the augmented generator. And finally, data were taken with the connecting 

<r~ 

arm shortened, to oscillate the flap unidirectionally from = 0° to 
^ f- = 60° for dynamic observations -discussed later. 
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DISCUSSION OF EXPERIMENTAL RESULTS 



STATIC : 



In the static case the experimental trends were correlated with the 

/■' f, 

theory . In Figure 8 a typical central position of J ^ was chosen and 
observed. Diffusion of the vortices downstream is negligible from a 
theoretical as well as experimental point of view. Westmore (Reference 29) 
relates a variation of with downstream distance as follows: 



= -a 





By substitution it is observed that, with the tunnel "q T s ,T in the present 
discussion and the short length of streamwise distance, diffusion down- 
stream is negligible. Even with the additional turbulence introduced by 
the tunnel walls, no dissipation can be observed over this distance. In 
fact, as is observable in Figure 7> there is a strengthening in ^ as the 
probe is translated downstream. This is most likely due to the decrease 

and increase in cos Ca for positions farther downstream when 1 r r 

(Figure 19). 



Having assumed that all of the vorticity is wrapped up in a single 
trailing vortex, the profile of dcwnwash (Reference 18) is as follows: 
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Although this graph is induced angle of attack with distance X, it gives 
a good indication of the actual variation with distance of the vorticity 
as a combination of circulation and jet sheet components to the total 
value. This, therefore, is a good representation of the strength of the 
trailing vortices. Here then we can see that as X increases it is obvious 
that at some central position J'^pj will be greater than at an even farther 
downstream. Observing the actual tendencies in Figure 7 we see 
that this hypothesis holds true in all cases. In the x = 2.0c position, 
the large value of 1\ and small cos Or, magnitude, due to core position, 
will overcome this initial ['at x = c+' hence, the value of uJp will peak 
at x = 6.0c downstream. This discussion represents the static case with 
the flap deflected positively, (toward the centerline of the tunnel). It 
also introduces the interference tendency far downstream due to wall 
effects that are not noticeable near the flap. 



When the flap is directed in a minus direction, (toward the closer 
tunnel wall), K will increase with x as will cos . Thus, their 
effects tend to cancel in the generation of top. Again P will vary as in 
Figure 6. Theoretically, the general trend for /‘Opj should be decreasing 
in value with downstream, distance. This is normally true (Figure 19). 



In the set of Figures 8 to l6, the entire cross-section is vi ex-red. 
Both magnitude and direction of the -^vectors make the position of the 
core in the static condition easily located. The solid vectors are for 
+ ^ f and the hatched vectors for - Sp . Magnitude and direction 
differences between the hatched and solid vectors illustrate the movement 
of the equivalent vortex core position across the tunnel. 



In any cross-section, movement obviously decreases 'Of, 

I 

with increasing \\ proving experimentally that 
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r 

Velocity and changes alone can be observed by comparing Figures 8, 
9 ; and 10. These are typical arrays of cross-sectional induced velocity 
vectors with varying free stream velocity . Thus, as is* evident from 
these figures, 




) 



The general trend is toward an increased for an increase in V 
however, the partial derivative value is not strong. Hence, this variation 
.is somewhat erratic. Experimentally, (Figures 17 and 18), this solution 
agrees with the theoretical evaluation in Appendix I. 



Flap deflection changes, have the most obvious affects on 

of all the parameters. As can be seen in Appendix I, 




All three add a positive increment to the partial derivative such that 





A comparison, in progression, of Figures 8, 11 and l4 will give a typical 
view of the change in magnitude of the induced velocity as a functipn of . 
Figure 20 gives a typical graphical picture of this progression. It can 
be observed that the seems to be of the order of an increasing 

power function. It appears then that in fact the partials 
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FIGURE 10. 
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FIGURE II. 
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FIGURE 15. 
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FIGURE 20. 
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DISCUSSION OF EXPERIMENTAL RESULTS 
DYNAMIC : 

The dynamic experimental data for the majority of cases agree with the 
theoretical trends discussed in Appendix I. In all regimes, the decrease in 
induced velocity amplitude with high frequency is evident. 

Induced angle of attack holding all the variables but constant^ was 

recorded as a function of reduced frequency (Strouhal Number), K. As can 

. y* 

be seen from Figures 22 and 23, magnitude peaks in tkp arise at intermediate 
values of K for those sectional positions investigated. As K is increased 
these data tend to show the same decreasing trend as Figure 21. 



The rapid decay of magnitude of induced velocity as a function of K in 

Figures 21 - 23 at high frequency shows strikingly the dynamic interference 

affect of high frequency oscillation of the airstream. Some of the high 

frequency data, however, does show a reduction trend about a biased value 
t K* 

of wp; this is due to the vortex filament movement downstream. 



In the augmented case, as is stated in the theory of Appendix I, at 

high K, the magnitude of , for the sinusoidal symmetrically-oscillating 

flap does level off to a constant value in the range Q i\Sr> , < } Aj" 

' /tJr** 3 % i < 



In Figures 21 - 23, a comparison was made of the induced flow 
generated without blowing augmentation. It is obvious that the trends of 
peaking with intermediate values of K and decay with high values of K 
occur here, too. But as before, the use of blowing creates a great increase 
in magnitude even ’at high K values. 






Experimental results (Figures 2^-27) indicate unidirectional flap 
oscillation decreases the cancellation of the harmonic components which 
appear in the symmetrical case. However * the theory proposes that this 
high frequency cancellation should be exactly the same"as that for the 
symmetrical deflection. The theory and experiment disagree at this point 
since the does not diminish symmetrically around the mean value. 

Vortex filament motion as described in Appendix I will only change 
the value of the bias such that at uj- = the j/J^Jwhen the filament is 
assumed motionless } is less than jc^pj when the filament is in motion. 

This fact is still of less importance than the previous disagreement 
between the calculated and the experimental cancellation. 

Although further experimentation at higher K is necessary to verify 
the experimental decay of ; op with K, it would seem impossible to create 
gust by means of an oscillating Jet flap and trailing vertex system for 
Strouhal Numbers much in excess of 2.5. 
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FIGURE 21. 
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STREAMWISE POSITION 6.0c, SECTIONAL POSITION 40"xl2' 
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FIGURE 23. 



STREAMWISE POSITION 2.0c, SECTIONAL POSITION 40” x 24' 

8 f =0-460° UNIDIRECTIONAL 
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FIGURE 24. 



STREAMWISE POSITION 2.0c, SECTIONAL POSITION 30"x30 
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FIGURE 25. 



STREAMWISE POSITION 6.0c, SECTIONAL POSITION 40"x24 

Sf =0— +60° UNIDIRECTIONAL 
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FIGURE 26. 



STREAMWISE POSITION 6.0c, SECTIONAL POSITION 30"x30 

Sx =0 — + 60° UNIDIRECTIONAL 
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FIGURE 27. 



NUMERICAL EXAMPLE 



To compare the proposed theory of Appendix I with the actual ex- 
perimental data 5 a typical set of parameters is chosen to be analyzed 
in the static regime. 

The following are assumed: 

(1) l/A-'lC ’ such that 

X = 6.0c downstream 
Y = 30" 

Z = 30" 

* 

(2) ^ ^ = 45° , static in the positive direction 

(3) \‘ c- ; = 100 ft/sec with a choked 

f\ 

. * . \u. = . 84 



From Table I: 

(1) C u = 3.49 

( 2 ) = 39 .^° . 



(3) Filament angle (90° -)( ) = 8.07° 



From the Biot-Savart Law 




P( uart <9.) 

v. _ — ' 

< V7 K 





) — > 
=z 



fioy)r*L 



2.1 t>Cutr*'jrecJiM Cf’o'^-StcT.oof.L. 

oF 7 H C "5 ORV Vi* O DA'JM Pc S n>oF) . 
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The value of h when calculated from Figures A and B yields 

h = 1.88 ft 



Also assuming this value of h, then 



“S//J B 




Thus, the cos 




= .923 



A 



( 2 ) 



Substituting these values into (l) and assuming (cB-behind the flap 
along the jet sheet is constant^the induced velocity at point P is 
as follows : 




( 1 ) 



or 




C L Vc 0 r- ( \ + GOS 0 ^ 

■^7 



and finally 







(3) 



This value represents a single vortex in free air. Now the experimental 
wind tunnel results can be seen in Figure 11 , where the magnitude of 
induced velocity is 



of 
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This correlation holds true within the range experimental error, for 
the great majority of surveyed points. This would lead one to believe that 
some parameters in the wind tunnel are possibly changing the wall effects 
and cancelling them. In Heyson's report (Reference 12) the affects of both 
open and closed wind tunnel tests sections are analyzed. The overall effects 
of these two are opposite in sign. It also makes reference to the possibility 
of cancelling any wall effects by introducing a slotted or porous tunnel 
section. Something in the tunnel, which is seemingly obscure at this time, 
is simulating this porous tunnel "effect. Further study is necessary to 
discover the exact problems. 



ko 



PROPOSED MODEL 



The ideal gust disturbance would be a two-dimensional gust. It would 
seem reasonable that such a disturbance could be achieved by the use of 
superposition of vortices. This report proposes just such a system. Although 
previous research has generated different methods such as the Texas A&M 
pulsing jet for longitudinal gusts (Reference 27) and Giannini Controls 
Corporation oscillating wall jet gust simulator (Reference 6), the 
oscillating airfoil method used at NASA Langley (Reference 7) for generating 
a sinusoidally oscillating airstream seems best for a two-dimensional gust 
generator. 

As Reference 7 indicates , many problems arose in generating high 
frequency gusts in an attempt to simulate high speed flight conditions. 
Appendix I demonstrates the inertial load problems that would indicate 
selection of a simple flap mechanism with blowing augmentation to produce 
similar, if not greater, gust magnitude capabilities. 

Purely vertical or horizontal gusts can thus be generated by a 
superposition model. By oscillating the jet augmented-flaps either in 
phase or out of phase by l80° , these gusts can be produced. The vortex 
array should be patterned to generate as near a constant cross-sectional 
gust as possible. Thus, the gust seen by the model being tested will be 
nearly constant in a spanwise direction. As will be seen later, with the 
flaps oscillating out of phase by l80° , the most uniform gust will remain 
on the centerline while varying in magnitude and direction. However, when 
oscillated in phase, the filament movement will cause an unwanted s-ftift in 
this most uniform gust position. The complex velocity potential is, 




(i) 



4i 



for a single vortex, where 

■? - X t H c 




Since /“ 



a = location of the vortex in the plane. 




by substitution of (2) into (3) we arrive at 




By comparing ( 5 ) with (6) we see that 




( 2 ) 



(3) 



(4) 

(5) 

( 6 ) 



(7) 



( 8 ) 



(9) 



( 10 ) 



h2 



Substituting we arrive at the radial velocity 




( 11 ) 



and the tangential velocity 




This equation illustrates the inverse proportionality between the radial 
distance and magnitude of the vortex-induced velocity. 

Since the proposed system has an array of four vortices its velocity 
potential is assumed to be the following, 




The array would appear as in Figure 28. 









-v, 
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Such an array would produce a purely vertical oscillating gust. 

In producing a gust perpendicular to the flap chord, the dynamic 
movement of the vortex filaments at positions farther downstream 
will all be out of phase. However, even in a high frequency mode 
the lateral velocity components will all cancel in the central area 
of the tunnel. The position of the most uniform spanwise gust velocity 
will not move as a function flap deflection* This presents the best use of 
oscillating jet flap position, frequency and phase of oscillation, to date. 
If the jet flaps were mounted on the side walls of the tunnel a pure 
horizontal gust could be generated. 

From the complex potential previously described can be derived the 
conjugate velocities in the x and y directions. 

i! =■ ix 

>■?- 

Where u = velocity in x direction 
v = velocity in y direction 



By differentiating^the following exists 
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In separating the real and imaginary parts such that 




ff 






Using these equations, the theoretical magnitude and direction of the 
induced velocity can be shown. The movement of the vortex filaments is 
involved in this analysis. The equivalent vortex positions at x = 6.0c 
were assumed. This provides a means of comparison for this case with 
the actual experimental value. 

The dimensions of Figure 28 will be assumed since the experimental 
data will be of that configuration. With the gust generation devices 
mounted to the floor and roof of the tunnel, the flaps are synchronized 
l80° out of phase with the bordering flaps as shown in Figure 28. This 
will provide as even a vertical gust as possible. Figure 29 illustrates 
theoretically the magnitude and distortion of '0 which must be tolerated 
when testing under such a configuration. This distortion, for a 2-foot 
span model, is approximately .40$ from mid-span to tip. 

Figure 3C S is the comparable experimental result due to super- 
position. Here again, spanwise distortion, in vertical w is important. 
It is approximately 40$ over the same 2-foot span model. 



Creating lateral gusts would seem at first, easily accomplished 
by changing the phase of the flap oscillations. Figure 32, is an 
illustration of this configuration. The generators in the array are 
all of the same phase. However, the filaments all moving in phase 
back and forth are ignored. 

In Figure 32, experimental superposition shows that the center of 
the array of filaments as appears in Figure 31, will move a certain 
distance to the left as the filament particles move downstream. Hence, 
at maximum static deflection, at k - 6.0c downstream, a very erratic 
pattern appears in the centerline position. This is intolerable when a 
relatively constant cross-sectional gust is required. In the production 
of these gusts, both lateral and vertical, combinations of parameters 
could possibly be found and programmed to induce a uniform gust flow 
over the entire span of the model. In addition, a half span model could 
be mounted to the tunnel floor with the previous vertical turbulence, 
thereby simulating a lateral gust. This , however , is left for future work. 



THEORETICAL 
DISPLACED FILAMENTS 

C L =2.33 
V eo = IOO ^/ sec 
S f =30° STATIC 
X = 6.0 c 




— 

I w I =40 ft/sec 



FIGURE 29. CENTER LINE GUST DISTRIBUTION ACROSS 
THE SECTION. 



EXPERIMENTAL 
DISPLACED FILAMENTS 
C L =2.33 
Vco = 100 ft/sec 
Sf = 30° STATIC 
X = 6.0c 






Iwl = 40 ft/sec 



FIGURE 30. 



THEORETICAL 



C L =2.33 
Vcq =100 ft/sec 
8f = ± 30° STATIC 
X = 6.0c 






Iwl =40 fi/sec 



FIGURE 31. 



EXPERIMENTAL 



C L =2.33 

705=100 ft /sec 

S f =±30° STATIC 
X =6.0c 




— 

Iwl = 20 ft/sec 



FIGURE 32. 



CONCLUSIONS 



There are some basic conclusions that can be drawn from the 
discussion and analysis of the experimental data and theoretical 
calculations as presented previously. " 



1) Since only the flap and not the entire airfoil is being 
oscillated, the inertial loads at high frequencies can 
be decreased considerably. 

2) By properly programming the vortex generators, the array 
of airfoils can produce a uniform cross-sectional gust 
in the vertical or lateral directions. 



3) In the dynamic case the flap was oscillated both bi- 
directionally and uni-directionally. Theoretically, the 
former should produce an oscillation in induced velocity 
symmetrically centered about zero. At very high frequencies 
the dynamic components tend to cancel themselves. Experimentally, 
this was substantiated. Theoretically, the latter should be 

centered about a finite value such that when ~ 

7 J 





St /97 tc. 



However, the experimental results in this series of tests did not 
fully agree. The experimental dynamic peak to peak value did not 
decrease symmetrically about the constant value. However, 
frequencies corresponding to Strouhal Numbers in excess of 1.5 
must be observed and more data taken to substantiate this 
deviation from the theory. 



4 ) Vortex filament movement is a factor in determining the value 
of the constant value of Oj p when Iaj ^ 

With the present apparatus, the steady state value for the bi- 
directional oscillation will be displaced from the zero 
position. This, however, does not influence the fact that the 
time varying oscillation is cancelled as if the core position 
did not vary. Core movement in the uni-directional oscillation 
again simply biases the steady state value, but again, had no 
effect on the dynamic cancellation. The experimental data 
confirms these theoretical observations. Hence, in the 
production of a very high frequency, sinusoidal, time varying 
gusts, the vortex generator is of little use. However, the 
apparatus can be operable for a low to medium frequency band- 
width while taking into account the rapid magnitude decay at 
the upper band frequencies. 

5) Attempts to account for the presence of the tunnel wall by 

image vortex systems indicate that little or no wall effects are 
necessarily accounted for in the calculation of In the 

numerical example this fact was shewn for an average test point 
where it can be seen that the test section actually simulates 

a porous section, and, theoretically, needs little or no tunnel 
wall correction. Further investigation is needed to discover the 
reasons for this trend. 
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APPENDIX I 



JET FLAP ANALYSIS 



Important aerodynamic parameters of the jet flap were fixed due to 
the problem itself. First, an airfoil of a symmetrical cross section 

c* 

was desirable to accommodate the changes in circulation throughout a - p 
to iOf range which would in turn produce symmetrical gusts. The MCA 0012 
section was chosen for convenience because templates with a nine inch 
chord were already in existence, thus simplifying construction. 



Once the chord was fixed, the flap size had to be determined. 
Preliminary two dimensional work, described in Appendix II, was done in 
the Princeton 2" X 36 " smoke tunnel to arrive at a possible C L for the 
3-D model. Knowlton, in Reference l4, discussed methods of determining 
for a 2-D airfoil section by measuring the streamline deflection 
of an upstream position. 



n > r 

Various tunnel dynamic pressures, flap deflections, and ^ were 
observed with varying flap chord ratio. The tests were conducted at 
30-50 ft/sec. Thus, the data were for the range of Reynold ! s Numbers 
oK\^< (c» <4.5/,^ and, hence, this regime corresponds very well 
with the 3-D model. 



VORTEX FILAMENT 



From 3-D elementary airfoil theory without blowing argumentation, 



L - C. ^ S 



or substituting for a constant chord model. 







( 1 - 1 ) 



(i-ia) 



From the simple circulation theory of lift and the Kutta Joukowski 



theorem, 




( 1 - 2 ) 



such that lift per unit span is equal to the product of air density, 
free stream velocity and circulation. 



By equating (la) and (2) we arrive at 




(1-3) 



Eliminating terms by dividing through we arrive at a value for 
circulation as a- function of 
Hence, from (3) 

f - C L \ !&> C 

~~ o 



\Ju>) and c . 



Since 




(1-4) 



(1-5) 



in that C is a function of angle of attack and flap deflection angle 
with no blowing augmentation, we realize that when the angle of attack 
is zero then lift coefficient depends only on flap deflection or 




at o 



(l-5a) 



For the model used in these experiments, the flapchord is nearly 
half the total chord. Let us assume that the change in angle of zero lift 
is approximated by half the ^ p, or the change in angle of zero lift is. 





for cy-O 



( 1 - 6 ) 



Since, 



for any non-augmented airfoil, 




is a constant, 
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then. 



P = ^ e 



(1-7) 



Thus, for an unaugmented, plain-flapped airfoil at c<[ = 0 

P ^ \l^ 

until the flap deflection reaches a value of twice the stalling angle 
of attack, or when 

The average induced angle of attack, with a fully attached flow 

<r" 

is dependent only on ^ . Thus, 

^-£2 /V. d'f 

Vc& 

In the case of the blown flap augmentation (Reference 19), the lift 
coefficient is a function of both j p. and \)jP 3 where 



( 1 - 8 ) 



(l-8a) 



( 1 - 9 ) 



r - a 






i. 



( fr- ) 



( 1 - 10 ) 



p S' 

due to the introduction of y(. into the analytic determination of 
in Reference 20. 



Hence, circulation for the blown flap is proportional to the product 

s\ 

of '.1(0, and C u , with the latter a function of both L 'p andC^ . Thus, for 
this case 



D 



( Of CJ * 1 / ^ 

I v / ^ '-My v ^ 



<3 



(i-n) 



Wow by definition 



o* = ' " '' ' 

I s 



(1-12) 



in that the product of the mass flow rate of jet air and its local slot 
velocity is nondimensionalized by the mass momentum over the wing. Thus, 
is in fact, a lift coefficient in addition to the or circulation 
lift coefficient. ■ is constant throughout the present experiments; 
hence, dependent only on 

As in Reference 25, the increment in lift is by reaction to the 
vertical component of the jet momentum and the increases in wing 
circulation lift from the introduction of asymmetry to the airfoil by 
the jet and the energy injected into the boundary layer to retard 
separation. Both increases in lift will be assumed to add a component 
of to the trailing vortex shed from the tip of the 3-D 'wing. 

JET SHEET 



From thin airfoil theory, the jet must first be considered in 
inviscid flow • Figure 33 gives an illustration of the problem , 




The high velocity air jet is issued from the. leading edge of the flap , 
leaving the flap surface at the jp angle . Both the flow inside the jet 
and surrounding the jet are irrctational, Total head pressure is greater 
in the jet than in the free stream. Figure 3^- illustrates a segment of 
the jet downstream. 
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Since the flow has a discontinuous velocity gradient, we can consider 
the jet surrounded by two opposite signed vortices each proportional 
to the strength of the local flow external to the jet. 



Again from Reference 25, the pressure at the boundaries is 
continuous, i.e. total head is constant across the boundaries. From 
Bernoulli T s equations, 



O 





y u JD i ) 


(1-13) 




4 


(l-i4) 



Since irrotational flow is concerned 




(1-15) 



By involving mean velocities in the discussion and rearranging algebraically 
we can arrive at the conclusion that 




J 



(1-16) 



6 o 



where 




( 1 - 17 ) 



Again this assumes that the Jet extends to infinity and thrust is 
totally recovered* „ 

The pressure differential described in(l-l6)is due to the difference 
in the elemental strength of the upper and lower boundary vortices. 

By integrating the two into one vortex by assuming 



we will arrive at a vortex sheet with a strength equal to the 
difference in strength of the previous two, located at the centerline 
of the jet. Thus, 



Following Spence T s progression in theory, the strength of the vortex 
jet sheet per unit span is 



/- 



- O 



'j 





( 1 - 20 ) 



But since 
approach oo as 





( 1 - 21 ) 
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and the pressure differential per unit span can be assumed from (l-l6) 
to be 

_ o2j_ 

R. 



(1-22) 



Now the total lift is the sum of circulation lift and momentum 



lift 



n 






(1-23) 



n 



Where I ^ is the corrected circulation lift for the blowing augmented flap, 
and 



Since 



OC 



0 = i ds 



X- '= ^ 

J 2R 7 



(1-24) 



(1-25) 



Then the following thin airfoil approximations are made for small J /= 



and 



1) 



2 ) 



the vortex distribution representing the jet is located from 
I ^ /• < 0 ° along the real axis ; 



x-Cx) - 



j ^ ~ 

where y (x) is deflection of the jet below the horizontal line; 



and Y c (. Yy " ^ O 0 Y L^-) 



as in Figure 35 



3) 



for 0 ^ X ^4 / 
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(l-25a) 



By following through with integration of(l-25a), taking the function to 
its limits, and solving for C^a total of the force components will 
give values of 

cKi 3 5<i f 1 .. ^ -l Cc 4 - Cyj 



'7^ 



c) < 5 ) 



F 



(1-26) 



and 

_ p7 1 -t- yi 

d oi 

With further refinements Spence ' s equations are displayed in McCormick 
(Reference 20) as follows: 

S 6c / 



(1-27) 



and 



d £c 



( 1 - 28 ) 



on (' i +' / t ' 1 2/1 j/i 



(1-29) 



Q/ 

These equations hold true for / [^ = 0, however , the graphical plot in 

Reference 20 was used for actual interpolation between the above two 
equations. 
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For 



total u 



the 



following is true: 







W 1 



F 




(1-30) 



Since the present flap chord is in the region of 

0< Cr /c <10 

A 

an interpolation was necessary, and a numerical ^ was obtained. 

Table I lists the numerical values for all experimental regimes calculated 

using McCormick 1 s interpolation. Also the actual value of 

n f 

is listed for varying values of ^ and d f • Numerically, here it can 

r • i ^ 

be seen that' -(or v*?) has little effect on the average downwash angle 
at infinity; but that it is nearly linearly dependent on Jr. 



Thus 



while 



C 

d 

'a Ck \ 



iO 

■40 



(1-31) 

(1-32) 



with a fairly constant positive value. This in turn shows that the 
position of the vortex filament downstream, as shown in Figure 36, is 
not a function of \U? but only of Also this leads to the conclusion 

that pv and cos Gf, are independent of 



For analytic purposes assume the value of 
be correct. 



gCu given in (29) to 



Thus, 





o, 



\- 










r 

■V f'vT/A 




(1-33) 



where 



TABLE I 



V^(f 50 75 100 

3-33 1.50 .84 




U) 

0 

0 


3.4 


2.62 


2.33 


45° 


5.1 


3.8 


3.49 


60° 


6.8 


5.07 


4.66 


30° 


25.7 


26.2 


26.3 


45° 


38.6 


37-9 


39-4 


O 
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VORTEX FILAMENT POSITION FOR SINGLE VORTEX 




FIGURE 36. 
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Since from (4) 




Then substituting ( 29 ) and (33) into (34) we get 

i 1 ■= sv-ft ^ 

IZ'-J 4 

where 
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From the lav; of Biot and Savart, the induced velocity anywhere in the 
system due to the circulation in a single trailing vortex is, 
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From Dwinnell (Reference 4), in the steady flap condition, the total 
induced flow at a point "P" due to the entire trailing vortex is 
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Expressing T and in terms of ©and d© 
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After integrating (38) the following results: 
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For an infinite vortex 
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and for a semi-infinite vortex 3 
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In the present experimentation the value of induced flow at a point P 
will always be in the region 
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for the static case. 



<; 60 



P 



/ 



lit u 
^ i / r\ 



(1-42) 



Since 



© 5 •?- 1 sc 



,0 



Then COS 0 



k 3 



J 



6 8 



Hence from ( 25 ) 
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Now assuming a constant spanwise downwash distribution in the work of 
Spreiter and Sacks (Reference 26), we know that the trailing vortex 
core will flow in the direction of the flap deflection with a velocity 
of approximately 



VcdRC = .41 tan C<\ 



Vc>3 



or as in Reference 18 
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Thus, 

^ v rW?s 4 

Then by substituting (l-29) and. (l-35a) into (1-46) we arrive at 




1 - 43 ) 



1 - 44 ) 



1 - 45 ) 



1 - 46 ) 



( 1 - 47 ) 



Transposing (1-47) we arrive at 
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we can readily see that the value of ~^~\/Zo 

is independent of velocity and is dependent only on . This theoretical 

analysis corroborates the previous hypothesis that ‘^has no affect on 

either \ or cos^. 

/ - 
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Thus, the partial of circulation with respect to ’.(p approaches the value 
given as follows: 
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For velocities decreasing from 100 ft/sec to 50 ft/sec 



C_3_ 
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Hence, the partial derivative is a positive number, or 




A typical plot is shown in Figure 38, for P as a function of 0> . The 
actual value, though positive, has a small value and hence the affect 
of \(pon f is very small. 



Previously, it has been proven that 
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This too is obvious from the plot of 



in Figure 38. 






as a function of 
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From ( I— 1+3 ) it was shown that 
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Assuming (l-53) is true and substituting (1-54) into (l-55) it is 
obvious that ' is a function of >*obut is influenced very weakly by 



a change in C-. Thus, 
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and is approaching a constant at high '.Os . 
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FIGURE 38. 
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C* • r - 1 

Next, the influence of jj- on Liy can be investigated. Again from (l-55) 



we know that 
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(1-56) 



Readily, we can see from that (1-35) that 
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constant 



Geometrically, the. value of i can be described as 
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from Figure 39 j 
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Now the cos\^ is directly related to k in that 
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Assuming small angles throughout 
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Hence from (I- 65 ) 
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Now* recalling from (1-46) that 
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We know that the derivative of Jy(with respect to y)^c is 
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C /- 

Thus, y-fhas no effect on iyi. Now recalling (l-29) it is known that 
with all other parameters remaining constant 

6 S r - ~ c °' J ^ r A 

Hence ? from chain rule differentiation, the results are, 
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From (l-69) 
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So by substituting (l-7l)- and (l-73) into (1-72) the result is obvious. 
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Hence, i s linearly proportional to 

From both (I- 67 ) and (l-68) 

J D , CC&' 4 •C'* ^1^ 



oJci 



Hence , 



3 Qf\ 






4 



p- 



(1-71) 



(1-72) 



(1-73) 



(1-74) 



(1-75) 



(1-76) 



Both values in ( 1 - 76 ) and (l-77) are not constant. 



Incorporating ( 1-57 ) 5 (1-76) and (l-77) into (1-43) we realize 



that 




( 1 - 78 ) 



This is contributed to by all three variables defined in ( 1 - 56 ). 



From the previous study we conclude that for the standard flapped 
airfoil without blowing augmentation that 6^'^ - j- whereas in the 

blown flap case it is (yJp> ~ ^ ^ 



In the regime studied, n<’l. At higher velocities it approaches 

f J V ^ Fee. /rv - / 

In summary it can be said that there are two variables of importance 
in the present study. They are j/r and All other parameters involved 

in the generation of induced velocities at remote points downstream are 
also functions of SV : and And, from the previous discussion it would 5 * * > 

that the more important parameter in the final magnitude of generated 

turbulence. 
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THEORETICAL DISCUSSION OF DYNAMIC INTERFERENCE 



As discussed in Reference 2U ? the forces on a flapped airfoil are 
a function of its frequency of flap deflection and its aerodynamic 
properties. The magnitude of induced velocities varies as in Figure 40, 
a vector diagram following the Theodors en function in magnitude and phase 
changes . 




As can he seen, the phase lags, and the magnitude decreases with increasing 
reduced frequency, K. In this discussion 



K = c 

\)o° 

In the experimental work the reduced frequencies vary from 



(1-79) 



0 ^ K < /. 5 



(1-80) 



Obviously, with such high values of K, the magnitude is decreased by 
nearly half, for the unidirectionally oscillated flap. 

Two unsteady effects have a significant bearing on the forces and 
moments on the wing. "Virtual mass" and "wake effect" are the component 
parts that vary the forces and moments on the wing/flap with increased K. 

The virtual mass effect has no influence on the strength of the shed 
vortex. It does, however, increase the load on the motor which oscillates 
the flap. By incorporating a motor with a large output, this problem is 
of no consequence. 
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The wake effect is a 3-D effect as opposed to the 2-D virtual 

mass problem. Although it can be shown mathematically, that the 

flow-induced by a flap oscillation through a plus and minus j would 

theoretically cancel itself out at very high the sinusoidal 

oscillation was attempted for a first run. The final value would be 

zero as opposed to the case shown in Figure 43 ? in which the flap is 

oscillated in one direction only; in this case^ (, Jp would have a constant 

value when [a) r # 

1 



For the symmetrically oscillating case the magnitude did begin to 
decrease when higher frequencies were encountered. However, it is 
obvious that the theory will be altered when the movement of the 
vortex core is taken into account such that at c/? , the 

actual was still greater than zero, although only a fraction of 
the static or quasi static value. Figure 4l, illustrates the wake 
effect for symmetrical oscillation. 



The law of Biot-Savart gives the incremental induced velocity 
at "P" in terms of circulation strength, ! 1 ; distance from the 
vortex core, h; and the difference in the cosines of the angles to 
the ends of the segment. At high frequencies, both positive and 
negatively signed vortex circulation components affect the induced 
velocity at P. Therefore, the actual magnitude in the positive 
direction is less than that of the static or low frequency case. 
Figure 42 shows the static case in which the 
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circulation is assumed constant to infinity and that there is no 
dynamic cancellation. 







By oscillating the flap in one direction, the only real affect on 
the sinusoidal frequency response is to input a bias such that as 
the value approaches zero. Hence, assuming at present that the 
vortex core position is constant, the oscillatory component of the 
induced velocity will symmetrically decrease about the mean value of 
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until at iO.-i d- 0 } (pj.will actually equal this constant value and the 
oscillatory gust will no longer be present. 



Figure 43 illustrates the typical core motion and accompanying 
circulation strength for a deflection in a single direction. 




Since the vortex core does cross the plane of symmetry it is obvious 



The sensor dynamics had little or no affect on the output data since 
the experimental oscillation bandwidt h of the present discussion was far 
below the natural frequency of the* sensor. As described in Appendix l¥, 

= .7 for the sensor* These frequency and damping characteristics 
assure fidelity in the measurement of the induced velocities. For 
experimentation at much higher frequencies (i.e.; 6^ >40 hz) sensors 
with correspondingly higher natural frequency must be substituted for 
the present vanes. 



that the bias will change and that the value at bj^^will be 
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APPENDIX II 



INITIAL EXPERIMENTAL WORK IN TOO DIMENSIONS 



Experimental trends in generation for the blown flap were 



arrived at through a comparison to an analytic and experimental study 
of a simple airfoil section with known aerodynamic properties. Some 
characteristics of the blown flap model had to be preserved. First, 
the slot position was investigated. From a physical point of view, the 
slot must be directed to produce an increase in boundary layer energy 
in the streamwise direction. This re-energization of the boundary 
layer will tend to maintain attached flow over the flap and 
thus from supercirculation addition, increase C L . 

Again, as previously analyzed, the vector component of momentum also 
adds to the . The slot is chosen at the flap hinge point since 
it is there that separation would be likely to occur and it must 
operate bi-directionally. Since y in all cases, the forward stagnation 
point will always occur exactly at the leading edge for 3(^0 



not necessary. Slot area was designed at a width capable of being 
choked at the required flow rate by the Schramm compressor. 

The chord of the airfoil and flap are important in lift and 
circulation development. To correlate the 2-D experiments with those 
of Hazen and Knowlton (References 10 and l4), the overall chord was 
chosen to be 16 inches over the tunnel velocity range, the Re T s 
agree well in both 2-D and 3-D cases in concern. The flap size 
is of great importance since the airfoil will always have an ^ 
initial 0 ; hence, the larger the deflected area, the larger the 
circulation produced. The flap size, however, must also be 
considered from an actuation load viewpoint . Thus, the required 
circulation must be weighed against the required input power to 
oscillate the flap at high frequency. A basic .45c flap seemed 
a reasonable initial compromise. 



and 




very large A * , leading edge blowing is 



The tradeoff of deflection angle for inertia of the flap 
was also a basic concern. Although very high flap deflection 
angles could produce high C L , it is also a fact that the blown 
air produces dynamic loads which increase with j f and add to 
the already large inertial loads at high frequencies. Previous 
experimentation by Korbacher (Reference 15) with jet flap blowing 
has lead to a group of optimal parameters for a similar model. 

In his work, Korbacher optimizes the lift and drag. In the 
present study, the drag is of little importance since it is 
circulation only that produces shed vortices. Thus, as large 
a deflection as possible -would be the optimum,. Actuation 
mechanics again limit the range of operations. 

Flap deflection greater than 60° would be difficult to oscillate 
with the present equipment; hence, was taken as the 

maximum with variable connecting arm positions available for ±45° 
and ±30° . 

2-D Smoke Tunnel Testing 



The 2-D smoke tunnel experimentation was used as a means of 
observing trends in production due to variations in flap chord, 
and . The plot on Figure 47 gives some idea as to the 

results of these tests. Knowlton in Reference l4 describes a method 

> > 

of calculation of L/_ by streamline deflection upstream. In the 
present study a simple comparison was made between the jet flapped 
models and the 2-D symmetric airfoil with assumed Q = .l/deg. 
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Having made this assumption, the 2-D symmetric airfoil without 
blowing was given increments of Oi . Observing the streamline 
deflection at a point K upstream 




of the leading edge in Figure 44. Thus a graph of versus -\< 
H I 

and hence, L L versus .^was arrived at for the perfect 2-D, 
infinite airfoil. 

The assumptions made to this point are: 



(l) infinite /Tv 




0 



(3) 




0 



(4) varying J>C 



(5) tunnel velocity ^ !o=>^ = 50 ft/sec 



The results of these observations appear on Figure 45. The 
versus Ko has a constant slope. Thus, this curve can be used 
to detennine ' for other configurations at C>(_ = 0 but with 





FIGURE 45. 
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Results of this comparative analysis appear in Figure 46 for various 
{J and S"V at = # 25c and then .35c. As can be seen, the overall 
trend illustrates that the various partial derivatives in ^production are 
all positive. From these conclusions and knowledge of dynamic inertial 
load difficulties 9 the following values were chosen: 

(1) C( /z 

( 2 ) - ± 

(3) CHOKED FLOW 



Vortex generation was attempted at NASA Langley (References 7 and 22) 
in the low speed wind tunnel facility. High frequency oscillation of 
the entire airfoil seemed mechanically unfeasible. By decreasing the 
size of the oscillating area to just the flap aided with blowing, the 
large C is possible without the tremendous inertia loads at high 
frequency- As discussed previously from Figure 46, the values of 0 



and 



' c 

^ t are obviously high. With c 9c increased to .45c and 

will increase greatly. 



to a choked condition the entire 



CONCLUSIONS OF 2-D SMOKE TUNNEL EXPERIMENTATION 



l) A basic symmetrical airfoil with a .45c flap will give an 
adequate tradeoff between ^ and the aerodynamic and 
inertial load effects at high . 



2) Choked flow at the blowing slot will give an even emission 
of high velocity air across the span of the model while 
producing the largest possible . 

- 1 60 will be the optimal tradeoff deflection between 
/v(\ and mechanical driving problems. 



3) 
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4 ) From the graphs (Figure 46 ) previously presented of 
experimental C_ production, partial derivatives can be 
assumed as follows: 





C 



4 



and 
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APPENDIX III 



/ 



SCALING 



In wind tunnel testing, scaling of parameters is very important 
in the simulation of real atmospheric situations. o£ major importance 
in producing acceptable turbulence in a 3-D tunnel is the frequency 
of the sinusoidal gusts being generated with respect to the tunnel speed. 
Reduced frequency, K •= Cd A is a nondimens ional parameter that 

expresses this relationship and represents a grouping of the important 
dimensional characteristics. 

To begin^we first hypothesize a full scale aircraft and a model 
size. If we assume a gust wave length equal to the overall length of 
the aircraft we would have typically 




and 






rroi> 




(III-D 



(III-2) 



Next the range of K possible for the wind tunnel can be calculated. 

— Off (/Xrf ? O n.co 



k' 

‘'rACD 
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V 



YtC i> 



Here 



X >.cA- ~ A Yc.C'ii -Vt ■ 



(m-4) 



5 0 < ioq 



with 



• (m-5) 
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t 

and 



O < kJ'hxCD ^ | ^ k-?-' 



(III-6) 



the possible range of K is 

0 < K < -4.01 



(III-7) 



Assuming 



and 



Va/c =600 Hfckc* 



4/c - 50 -ft. 



and since 






«/c 



/ c ' 



(III-8) 



(m-9) 



(III-10) 



for a correct simulation , 

w V - K nca {l«i) 

or . / 

0 < w fl/c <'76? A? 

The high frequency simulation is of interest for the study of 
structural and aeroelastic problems ? whereas the low frequencies may be 
of interest in simulating pilot-control inter-relationships in the 
stability and control region of aerodynamic study. 



hz (III-11) 



(III-12) 
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APPENDIX IV 



VANE DYNAMICS 







The dynamics of the vanes are very important to accurate 
measurements in both magnitude and frequency. Graham (Reference 8) 
undertook to hypothesize the dynamics of just such a vane with 
optimal dynamic characteristics. 

Vanes , connecting arm, and potentiometers are the essential parts 
to the angle measuring system. All contribute to the overall moment'' 
of inertia of the system and thereby affect the natural frequency 
and damping in the simulated second order system. 

Such a system can be described with the basic equation of motion, 



This equation describes the summation of torques due to angular 
acceleration, angular velocity and angular position. * 

First, the torque tending to bring the vane to zero angle of attack 
is 




(IV-1) 




This equation is the summation of the products of lift and moment arms 
for both fore and aft lifting surfaces. 



90 



For small angles 



Cod 0^1 



Hence, 



7 



ft 




(lV-2a) 



(iv-3) 



The surface area of both vanes is identical so 

[i» -*f~] 



Now 



o = 9 j><< c< ^ 





In both cases 



o(=Q 



Thereby 
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(I V-6) 



(lV-6a) 



(rv-7) 

(iv-8) 



Substituting in the torque equation 
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(I v-9) 
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Differentiating with respect to 0 yields 



(p- - ) C-* KdoJ^ ^ ^ 

d.a ^ 1 0 



(rv-io) 



For small angles of attack 
and C_t) c I s approximated by 



& VT\ I tsj 



Thus 




(IV- 11) 



Next, the torques proportional to angular velocity are important. 
As angular velocity occurs a damping opposing the motion is 
developed from the change in velocity perpendicular to the free 
stream. 
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Differentiating with respect to O we see 
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(IV-12) 



(IV-13) 



Collecting terms yields 



Jie 




n - Q ^ 

■V^ 




(iv-i4) 



Substituting (lV-l 4 ) and (IV- 11 ) into (iV-l) we arrive at the 
characteristic equation of motion for this second order, mass-spring- 
damper type system. 




4 - 
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de 




or 





(iv-15) 




- ,0 



(IV- 16) 



By comparing coefficients of (iVt 16 ) and (IV-l) we see that 






i] 



(IV- 17) 



and 



y 



u = \ 



J 



C^ACr>„SJ <&-'?) 

\L 1 



For the experimental model being used and the dimensions of the vanes, 
the following dimensional characteristics apply: 

2= dU, V 

„ G 0 ( 
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Solving equation (IV-17) 

uv =55 Jk 



and 







The following vane system is therefore critically damped with a natural 
frequency of 35 hz as indicated. 

Experimental tests were made and the vanes react exactly as the 
theory would predict. Figure 49 is an example ^of the vane reaction 
to a step input. 

It can be concluded that since the of the vanes is far above 

that of the actual testing frequencies and since the vanes are 
critically damped, no overshoot will occur even at high frequencies 
and measurement accuracy in the dynamic case is commensurate with 
that in the static case. 
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